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1. INTRODUCTION 

Advanced mobile systems aim on supporting high data rate services over multiple systems for multi- 
media applications involving voice, data, pictures and video. They should provide services without 
interruptions for private and public networks. These services require interference free transmission over 
frequency selective fading channels. OFDM is a promising technique for Long Term Evaluation-Advanced 
(LTE-A), Worldwide Interoperability for Microwave Access (Wi-MAX) systems, IEEE 802.11ac, and IEEE 
802.15 [1]. A good channel estimation procedure is necessary in order to make reliable decisions and to 
coherently demodulate the received data. OFDM systems require effective and exact channel estimation for 
time-varying channels. 

The principal motto of this paper is to improve the OFDM system performance by proper channel 
coding, suitable digital modulation schemes, adaptive filter channel estimation techniques and multiple 
antennas at the transmitter and receiver. By using all these, an integrated OFDM system is proposed to 
enhance system performance [2-3]. Gosh [4] formulated a closed form expression for impulse noise 
performance in a multicarrier QAM system and showed that it is better than single carrier systems for 
moderate values of impulse power and probability. Dong et al. [5] determined the power and rate adaptation 
analytically with channel state information in the presence of carrier frequency offset (CFO) for an OFDM 
system. Sabbaghian et al. [6] developed a new code based on dual BCH code for reduction of PAPR in multi 
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carrier system. Hong yu et al. [7] developed an iterative based ML algorithm for multiple solution 
complication from log likelihood function in OFDM mobile wireless channels. Kalyani et al. [8] proposed a 
channel frequency response estimation algorithm with limited number of pilots to determine the selectivity of 
the channel in an OFDM system. Han oh et al. [9] described an estimation scheme to detect symbols which 
are determined by the adaptive threshold using the noise variance with the different posteriori probability at 
each sub-carrier. Gao et al. [10] developed adaptive frequency domain channel estimation for uplink in 
MIMO-OFDM system to improve channel estimation by making use of correlation between two sub-carriers. 
Souryal et al. [11] reported that the sub-channel bit and power can be allocated dynamically to the state of the 
channel rather than fixed and uniform modulation. It was shown that bit allocation based on the prediction 
state of the channel decreases the BER by a factor of two with a maximum Doppler frequency of 50 — 100 
Hz. The performance of adaptive OFDM over non adaptive OFDM is good, even in the presence of channel 
estimation errors for slowly varying channels. 

Hong et al. [12] proposed adaptive all pass filter scheme in which phase rotation is used to reduce 
PAPR in MIMO-OFDM systems. Wang et al. [13] modeled adaptive receiver that is used for reduction of 
inter carrier interference with the help of optimal phase rotation. Elangoven et al. [14] described a pilot 
channel estimation in which fading effects in the channel are equalized by using LMS, NLMS and RLS 
algorithms. A Salh et al. [15] described the achievement of higher data rates through zero forcing beam 
forming with less number of radio frequency chains. Erik G. Larsson et al. [16] stated that the spatial 
multiplexing is used in Massive MIMO to increase the capacity, nearly ten times or more and improve the 
radiated energy in the order of hundred times simultaneously. Bolcskei et al. [17] described the fundamentals 
of MIMO-OFDM wireless systems and hardware implementation aspects of fourth generation cellular 
systems. Stuber et al. [18] discussed issues related to synchronization and channel estimation in MIMO- 
OFDM along with space coding techniques. In [19], serial concatenated convolution code for MIMO systems 


are described under Rayleigh fading channels. 


2. SYSTEM MODEL 
Let N be the number of subcarriers in an OFDM symbol, then the baseband discrete time OFDM 
signal sequence can be written as 


j2mkn 


x(n) = IFFT(X(K)) = o X(K)e N (1) 
k=0 


Where, n= 0,1,2,...N-1 and K is the index of the sub-carrier. Before transmission, cyclic prefix (CP) is 
appended periodically to prevent x¿(n) from ISI. The resultant OFDM symbol is transmitted through 
frequency selective fading channel h(n). Noise is added to the receiving signal in the channel. The received 
signal is expressed as 


Yg (N) = xg (N) * h(n) + v(n) (2) 


Where v(n) is the additive white Gaussian noise. At receiver, serial data is converted into parallel and CP is 
removed from y,(n) to form y(n) and is expressed as 


N-1 
jamnk 


Y(K) = FFT} = X y(ne 3) 


where k= 0,1,2,...N-1. Further, Y(K) can be written as 
Y(K) = X(K) H(K) + VK) k=0,1,..N-1 (4) 


Pilots are extracted from the obtained sequence Y(K). The estimated output X,(K) is obtained by: 


A AK) =a a (5) 


Transmitted data is recovered at the receiver by using proper demodulator. 
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Table 1. Parameters used in the Proposed System 
Parameter Specification Parameter Specification 
Number of sub- 


: 64 Coding Convolution 
carriers 
FFT size 64 Code rate 1/2, 1/3 
Guard type Cyclic extension Wireless channels Rayleigh, Rician and AWGN 
Cyclic prefix 16 (1/4) Channel estimation LMS, RLS 


BPSK, QPSK, 16- 


Modulation schemes QAM, 64-QAM 


Antenna configurations SISO, 2x1 MISO, 2x2 MIMO 


2.1. Characteristics of wireless channel 

In real time applications, the data transmission is much more complex due to many factors such as 
line of sight communication between transmitter and receiver, user’s density, multipath propagation, Doppler 
shift and so on and so forth. Hence, different channel models such as additive white Gaussian noise, Rayleigh 
and Rician channels are formulated. AWGN channel model is the most commonly used channels in 
communication systems. The density function for Gaussian distribution is represented in [20] as 


1 Gap? 
a ee -=z e< x< a (6) 

However, in this channel, multipath fading, frequency selectivity and interference are not considered 
and moreover, it is not suitable for terrestrial applications. Rayleigh and Rician channel models are 
developed for multipath propagation applications. The received signal is distorted due to attenuation, 
reflections, refraction, and diffraction during transmission. The statistical time varying nature of received 
signal power is described using Rayleigh and Rician fading distributions. The envelope of received flat faded 
OFDM signal in non-line of sight (NLOS) scenario resembles Rayleigh fading distribution whose probability 
density function (PDF) is characterized in [20] by 


2 


x x 
po) = ew (-7a)o> 0x20 (7) 


When a strong signal component exists in addition to multiple weak reflected components, the received 
signal envelope is in the form of Rician distribution [16] and its probability distribution function is expressed 
in [20] as 


ox (x? + A?) Ax 
EE ae (- 20? Jr (2) (8) 
forx = 0;A> 0 


Where I is the modified Bessel function of the first kind and zero" order and p is the mean due to fixed 
scatter or LOS path. K-factor of Rician distribution is the square of the peak amplitude of the dominant signal 
to twice the variance. K-factor is typically in between 1 and 10 and Rician distribution approaches Rayleigh 
whenever K tends to zero 


2.2. Channel coding 

The aim of the channel coding is to enhance the reliability of the communication system through 
improving bit error rate (BER) performance by appending redundant bits to the transmitted data. These 
additional bits are used to make the system capable of error detection and correction. In this work, 
convolution code with code rates 1/2 and 1/3 is considered for encoding the data. Convolution coding has 
memory that uses previous bits to encode or decode the input data stream. Figure 1, the convolution encoder 
has three variables: number of input bits ‘k’, output bits ‘n’ and memory block length ‘L’. The received bits 
can be decoded using a soft input-hard output Viterbi decoding algorithm. Soft input decision action is taken 
based on whole input range values while hard output decision action gives decoded bits based on fixed set of 
values. Viterbi algorithm uses code trellis and minimum hamming distance. 


2.3. Channel estimation techniques 

The subcarriers in OFDM system exhibit their own characteristics during the transmission. Hence, 
adaptive channel estimation is suitable for proper recovery of the transmitted signal at the receiver. Channel 
properties can be estimated using pilot based approaches. Block and comb type pilot arrangements are 
suitable for slow and fast fading channels respectively. In block type pilot arrangement, pilots are inserted in 
specific OFDM blocks periodically whereas in comb type, pilots are inserted uniformly in every OFDM 
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block [21]. Popular channel estimation schemes are Least Squares (LS), Minimum Mean Square Error 
(MMSE), Least Mean Squares (LMS), and Recursive Least Squares (RLS). Pilot arrangements for channel 
estimation as shown in Figure 2. 


4 - 4 
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Time 
Block type pilot arrangement Comb type pilot arrangement 
Figure 1. Convolution encoder with code rate 1/2. Figure 2. Pilot arrangements for channel estimation 
2.3.1. Linear channel estimation schemes: 
The received signal can be expressed as 
Y = XH +V (9) 


where X and Y are the transmitted and received signal vectors respectively. Here, H is the channel response 
while V is the noise vector. The LS estimator for an OFDM system can be written as 


ll 


Hs = (ex ay (10) 
is (11) 


The estimate of the channel at m" pilot subcarrier is given by 


Y,(m) 
H (m) = = (12) 
E Xp (m) 
The mean square error from the received signal can be described as 
J(e) = E[(H — #)?] 
(13) 


J(e) = E[(H - A)#(H - A) 


where ‘A’ is the MMSE channel estimate and X" is the Hermitian matrix of X. MMSE equaliser sets 
orthogonal by estimating vector Y as: 


E|[(H —H)y"]=0 (14) 

E[HY"]-—E[yy"]=0 

Ray — Ryy = 0 (15) 
The time domain MMSE estimate of ‘h’ is represented as 

f = RnyRypY (16) 


2.3.2. Adaptive filter channel estimation schemes: 

The LMS algorithm is a gradient-based procedure of steepest-descent that adapts the iterative 
procedure in which successive weight vectors are modified to obtain MMSE. The LMS estimator uses 
adaptive filter as shown in Figure 3. The initial weight is assumed and the successive values are computed 
based on previous estimation. The LMS algorithm iterates over each coefficient in the direction of 
approximate gradient. The weight vector equation is expressed as: 


w(n + 1) = w(n) +u [e(n)x(n)] (17) 
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where p is the step size, L is the number of iterations, e(n) is the error between received signal y(n) and 
reference signal d(n) . Here, x(n) is the input signal to the LMS filter. The minimum mean square error is of 
the form 


e*(n) = [d(n) - y(n)]? (18) 
The corrected received signal at the output of LMS filter can be shown as 
y(n) = w(n) * x(n) (19) 


The LMS estimate has lower computational complexity. However, it has poor tracking ability with slower 
convergence rate. In the RLS algorithm, previous samples are used to estimate the inverse of autocorrelation 
of input vector. The RLS algorithm uses past input data samples to estimate the inverse of the autocorrelation 
matrix of the input vector. A weighting factor for each sample is used to decrease the influence of very old 
input samples. The RLS estimator uses adaptive filter as shown in Figure 4. 


LMS efn) ain) RLS Adaptive | 4(") e(n) 
Adaptive filter filter + 
y(n) y(n) 
Figure 3. LMS Adaptive Filter Figure 4. RLS Adaptive Filter 


The RLS algorithm needs auto-correlation and cross-correlation matrix of input and is significantly more 
computationally complex while updating each coefficient, since it needs inverse of the correlation matrix 
¢(n). However, the RLS channel estimation algorithm has faster convergence with excellent tracking ability 
compared to LMS algorithm. RLS algorithm is described as; 


P(n) = 9" \(n) (20) 
E A *P(n — 1)x(n) 

k(n) oe en 

e(n) = d(n) —y() (22) 

w(n+ 1) = w(n) + k(nje(n) (23 


Here, A is the forgetting factor and for 0 < à < 1. The influence of past error values decreases exponentially 
and it is termed as ‘method of exponentially weighted least squares’ [14]. 


2.3.3. MIMO-OFDM 

OFDM transforms a wideband frequency selective channel into a number of flat fading channels 
which is suitable for MIMO structure that requires a frequency non-selective characteristic at each channel 
when the transmission rate is high. The MIMO is in the form of MrxMr, where My is the number of 
transmitting antennas and Mp is the number of receiving antennas. Let M;=2 and Mp=2, then the MIMO 
channel matrix ‘H’ can be expressed as 


h h 
H-| 11 d 24 
ho, hz C4) 


Here, ‘hj’ is the fading coefficient that indicates the impulse response between i"transmit antenna and j" 
receiving antenna. The MIMO system model is modeled as 


Y=HX (25) 


fled eee ca Ped (26) 
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MIMO is known as ‘spatial multiplexing’ as the received signal samples through all possible paths 
from multiple transmitting antennas are combined by the spatial ratio combiner at the receiver. This multiple 
transmit-receive antenna pair offers a linear increase in transmission rate for the same bandwidth without any 
additional power required. It can also increase diversity gain by transmitting the information signal through 
multiple antennas using space time coding. Alamouti scheme is the basic space time code for MIMO system 
is described in (27) 


= E | (27) 


where S, is the complex conjugate of S4 symbol and S," is the complex conjugate of S symbol. The encoder 
outputs are transmitted in two consecutive time slots from two transmitting antennas. In first time slot, S; and 
S are transmitted and in the consecutive time slot -S, and S, are transmitted. Received signals at first and 
second time slots are y(1) and y(2) respectively can be written as 


yA) = [ey hal [st] +n) (28) 


—S3 29 
y= [hy hal | 5] +2) e 
where h; and h, are channel responses in the 1“ and om path respectively, and n(1) and n(2) are additive 
white Gaussian noises at first and second time slots respectively. Hence, the MIMO system model can be 
described: 


be) ~ ns Sele (30) 


The column channel matrices C, and C, are represented as: 
hy hy 
C, = hl andC, = k 
The inner product of two columns C; and C, becomes zero as shown in equation (31). 


CEC, = hih, + h,(—-hj) = 0 (31) 


In equation (31) confirms the orthogonality of Alamouti scheme. The capacity and error rate analysis of 
MIMO communication system for satellite applications are well described in [22]. The capacity of MIMO 
channel is given in (32) 


SNR 
C = log>det (tv, + HH") (32) 
T 


where, I, is Ngx Np identity matrix, H is the channel matrix, and Nis the number of transmitting antennas. 


3. PROPOSED METHOD 

The proposed MIMO OFDM system is shown in Figure 5. In the proposed system, information data 
are encoded using a convolution encoder with code rates 1/2 and 1/3. The encoded data are applied to digital 
modulator in order to convert data into symbols and then pilot symbols are inserted for channel estimation. 
Then IFFT converts these symbols into time domain and makes the symbols orthogonal to one another. 
Cyclic prefix (CP) is appended to the output of IFFT to prevent the signal from ISI. In general, the cyclic 
prefix is taken from the tail of the previous OFDM symbol. The CP length should be greater than the channel 
response in order to remove ISI completely. After the insertion of CP, the OFDM signal is applied to the 
space time block code (STBC) encoder for transmission of signal via multiple antennas. CP is removed at the 
receiver and channel is estimated using either LMS or RLS algorithm. The reverse process of the transmitter 
is applied in order to recover the transmitted information. 
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Figure 5. Integrated MIMO OFDM system 


4. RESULTS AND ANALYSIS 

The proposed MIMO-OFDM system is implemented with the parameters mentioned in Table 1. In 
this paper, the OFDM system performance is analyzed for various channels as shown in Figure 6 and it 
illustrates the better BER performance of additive white Gaussian channel in comparison with the Rayleigh 
and Rician channels. Figure 7 shows the performance of OFDM system for various modulation schemes 
under Rician channel and it describes the performance improvement with decrease in order of modulation. 
Figure 7 and Figure 8 shows approximately 2 dB SNR improvement in LMS algorithm over RLS algorithm 
at 10-3BER under BPSK modulation. Better performance of BPSK modulation over QPSK modulation can 
be observed from Figure 9 to Figure10. The performance of convolution code with different code rates is 
analyzed in Figure 11 and Figure 12 which illustrates the improvement of system performance with decrease 
in code rate. Performance improvement of MIMO can be observed from Figure 12 and Figure 13. Hence, the 
proposed 2x2 MIMO OFDM system shows superior performance as depicted in Figure 14 over other 
combinations. Moreover, the proposed system performance is better over the system [23] in which the 
performance of MIMO-OFDM in terms of bit error rate is discussed for LS and MMSE channel estimation 
schemes. In addition, this system results superior performance over [24] in which the BER performance of 
MIMO-OFDM system for NLMS and DD-EW RLS channel estimation schemes is analyzed. Performance 
comparison of 2x2 MIMO OFDM system for digital modulation schemes under Rayleigh and Rician 
channels as shown in Figure 15. Figure 16 shows the system capacity improvement with order of MIMO. 
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Figure 6. Performance of 16-QAM OFDM under Figure 7. Performance of BPSK OFDM under 
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Figure 8. BER performance comparison of Figure 9. BER performance comparison of 
convolution coded OFDM using RLS convolution coded OFDM using LMS 
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Figure 12. OFDM system performance comparison 
for different modulation schemes and code rates 
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Figure 11. Performance of OFDM with 1/2 rate 
convolution code, LMS, RLS, QPSK modulation 
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Figure 13. Performance comparison of SISO, 2x1 
and 2x2 MIMO OFDM systems. 


BPSK 
QPSK 


BER 


~~ Ricial 
*— Ricial 


— Rayle 
— Rayle 
— Rayle 


—— Raylei 


160AN 
MQM 
gh BPSK 
hP | 
igh 160AM 
igh QAM 


SAR (db) 


Figure 15. Performance comparison of 2x2 MIMO 
OFDM system for digital modulation schemes under 


Rayleigh and Rician channels. 


Int J Elec & Comp Eng, Vol. 8, No. 5, October 2018 : 3829 — 3838 


Int J Elec & Comp Eng ISSN: 2088-8708 O 3837 


Capacity Vs SNR for different MIMO configurations 
T T T 


—-— SISO(1X) 
—S— SIMO(2X1) 
—— MISO(1x2) 
—*— mMIMO(2x2) 


a a a a a anaa aaaeeeaa 


15 


Capacity (b/s/Hz) 


o 5 10 15 20 25 
SNR(GB) 


Figure 16. Capacity improvement with increase in MIMO configuration for MIMO OFDM system 


5. CONCLUSION 

In the proposed work, the MIMO-OFDM system is simulated with 64 sub carriers, digital 
modulation techniques (BPSK, QPSK, 16-QAM and 64-QAM), convolution code with code rates 
(1/2 and 1/3), channels (AWGN, Rayleigh and Rician), adaptive filter channel estimation schemes (LMS and 
RLS) and MIMO configurations (SISO, 2x1 MIMO and 2x2 MIMO). From the results, it is observed that the 
BER performance of AWGN channel is superior compared to Rayleigh and Rician and also spectral 
efficiency increases with order of modulation with marginal degradation in performance. There is an 
improvement in system performance with decrease in code rate and the LMS channel estimation is found to 
be better than the RLS channel estimation. By combining channel estimation and coding, there is further 
improvement in OFDM system performance over individual schemes. 2x2 MIMO OFDM system shows 
improved performance over SISO and 2x1 MISO. Moreover, an increase in the order of MIMO configuration 
increases capacity of the MIMO system. Finally, the proposed channel estimated 2x2 MIMO OFDM system 
shows significant improvement in BER performance over other combinations. 
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